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ABSTRACT 

PSR BOS 18— 41 is one of the few pulsars which show multiple drift regions having well 
defined phase relationship. In this paper we report new results from the multifrequency obser- 
vations of this pulsar with the Giant Metrewave Radio Telescope (GMRT). We determine the 
mean flux of the pulsar at multiple frequencies and find evidence for a low frequency turn- 
over in its spectrum. Significant linear polarization is observed with depolarization at the edge 
of the profile, which is likely to be due to orthogonal polarization mode jumps. Circular polar- 
ization changes sign near the middle of the pulse profile at 1060 MHz which is not observed 
at 325 and 610 MHz. Remarkable frequency evolution of polarization angle (PA) is observed 
for this pulsar, which is commonly not seen in other pulsars. Based on the frequency evolu- 
tion of average profile, observed PA swing and results from subpulse drifting, we propose two 
possible emission geometries, a ^11°, (3 ^ — 5.4°and a ~175.4°, /3 ~ —6.9°. Simulation 
of the pulsar radiation pattern with both these geometries reproduces the observed features in 
the drift pattern quite well. In addition to the remarkable subpulse drifting observed at 325 
MHz, we report subpulse drifting at 244 and 610 MHz. At 244 MHz we observe drifting in 
the outer regions but not in the inner region whereas at 610 MHz subpulse drift is observed in 
both inner and outer regions. Sometimes, we see changes of drift rates, transitions from neg- 
ative to stationary and stationary to negative drift rates, many of which appear to have some 
connection with nulls. We determine (the measured time interval between the recurrence 
of successive drift bands at a given pulse longitude), P™ (the measured longitude separation 
between two adjacent drift bands) and A<1> (subpulse width) values at different frequencies for 
both inner and outer drift regions. Though P™ is the same for all the drift regions, P™ ™d 
A$ values are different. The peak emission from the leading and the trailing outer regions 
are offset by 9 Pi . Utilising this information, we solve the aliasing problem and argue that 
drifting is first order aliased with corresponding carousel rotation period P4 ~ 10 s, which 
makes it the fastest known carousel. The drift pattern from the two rings are phase locked for 
PSR BOS 18— 41. The same is found to be true for, all pulsars showing drifting in multiple 
rings of emission. This finding puts constraints on the theoretical models of pulsar emission 
mechanism and, favors a pan magnetospeheric radiation mechanism. 
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1 INTRODUCTION 

Pulsars with drifting subpulses are considered as an important key 
for unlocking the mystery of how radio pulsars work. Though 
there has been significant progress both in high quality observa- 
tions of drifting subpulses and in confronting their observational 
properties w i th exi s ting models (e.g.iD eshpande & Rankin ( 19990, 
iGupta et al.1 ( l2004l) . IWeltevrede et aU ( l2006l) . lBhattacharvva et al. I 
( l2007h (hereafter Paper 1)), the detailed nature of radio emission 
processes of pulsars and the exact location, distribution of the pulse 
emitting regions are still shrouded with mystery. 



Knowledge of the viewing geometry is essential to understand 
the pulsar emission mechanism. The orientation of the pulsar ra- 
diation beam can be determined by two angles, a (the angle be- 
tween the rotation axis and magnetic axis) and /? (the minimum 
angle between the magnetic axis and the observer's line of sight 
(LOS)). Possible geometry of a pulsar is determined by the "Rotat- 
ing Vector Model (RVM)" jRadhakrishnan & Cook] Il969l) ) fit to 
the polarization angle (PA) traverse. The possible choice of geome- 
try should also satisfy the average profile evolution with frequency 
and the observed drift pattern. 

The phenomenon of subpulse drifting, first reported by 
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brake & Craft] ( Il968l) . is manifested as systematic subpulse be- 
haviour: subpulses appear at progressively changing longitude in 
the pulse window following a particular path (known as drift 
band) which is specific to the individual pulsar concerned. Stud- 
ies of wide profile drifting pulsars can probe th e emis sion re- 



Table 1. Details of the observations at different epochs and conespond- 
ing measurements of the values of mean fluxes 



gi ons in the pulsar magnetosphere Gupta et al. Paper 
1, iBhattacharvva. Gupta & Gill ^20081) '). Since our LOS samples 



a large region of the polar cap in case of wide profile pulsars, 
presence of simultaneous multiple drift regions is quite probable. 
Phase relation between multiple drift regions are known o nly in 
few pulsars - PSR BOS 15+09 llMcLaughlin et al. I [iooi) . PSR 
B 1839-04 jWeltevrede et al. II2006I) and PSR B0818-41 (Paper 
1). Observed phase relationship between the drift regions can 
probe the distribution and dynamics of the emission regions in 
the polar cap. Multifrequency observations may offer insights 
into various aspects of frequency d ependence of pulsar radiation 
( IBhattacharvva. Gupta&Gill l |2008}) and the references therein). 

PSR BOS 18— 41 is a relatively less studied wide profile pulsar. 
In Paper 1, we reported the discovery of a remarkable drift pattern 
at 325 MHz, characterized by simultaneous occurrences of three 
drift regions: an inner region with flatter apparent drift rate flanked 
on each side by a region of steeper apparent drift rate. The drift re- 
gions have different values for P™ (the measured longitude separa- 
tion between two adjacent drift bands), but identical values for 
(the measured time interval between the recurrence of successive 
drift bands at a given pulse longitude). The closely spaced inner 
and outer drift regions always maintain a phase locked relationship 
(hereafter PLR). The observed drift pattern was interpreted as be- 
ing created by the intersection of our LOS with two conal rings on 
the polar cap of a fairly aligned rotator (q ~11°), with an inner 
LOS geometry (/3 ~ —5.4°). Results from simulations reproduced 
the observe d features (width o f average profile, single pulse drift 
pattern etc). IOiao et al. 1( ll995l) have studied the average polariza- 
tion properties of this pulsar at 660 and 1440 MHz. Their study 
reports a steep-sided wide profile with modest amount of linear po- 
larization (32% at 660 MHz an d 31% at 1440 MHz). From the wide 
integrated profile, [oiao et al. 1 ll995i) also predicted that this pulsar 
must be a fairly aligned rotator. 

We made multi epoch observations of PSR B08 18-41 at 157, 
244, 325, 610 and 1060 MHz using the GMRT in full polar mode. 
In Sect. |2|we describe the observations. We determine pulsar flux 
at individual frequencies and investigate the spectrum in Sect.|3] In 
Sect. |4| we illustrate the average profile evolution with frequency 
and in Sect. [5] we present a study of polarization properties of this 
pulsar. Possible choices of viewing geometries for PSR B0818— 41 
are discussed and compared with results from simulations of the ra- 
diation pattern in Sect.|6| Multi-frequency study of subpulse drift- 
ing is detailed in Sect.jT] Unique phase relationships between mul- 
tiple drift regions observed for PSR B0818— 41 are investigated in 
Sect.jSj In Sect.|9| we present a final summary of different aspects 
of the analysis results and the interpretations. 



2 OBSERVATIONS AND PRELIMINARY ANALYSIS 

The GMRT is a multi element aperture synthesis telescope consist- 
ing of 30 antennas, each of 45 m diameter, spread over a region 
of 25 km diameter. The GMRT can also be used in array mode 
by adding the signals from the individual dishes, either coherently 
or incoherently. We performed polarimetric observations of PSR 
B0818-41 at multiple frequencies (157, 244, 325, 6 10 and 1060 
MHz), with the GMRT, using the phased array mode jGupta et al. I 
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a : relative calibration 

b : absolute calibration 

t : simultaneous dual frequency observations 



iooS)), at several epochs. For t he polarimetric observat i ons w e 
used the technique described by iMitra. Rankin & Gupta ] (l2007h . 
The simultaneous dual frequency observations reported in this pa- 
per were done using the "band ma sking" technique described in 
IBhattacharvva. Gupta & GilU2008h . The signals from different an- 
teimas in the observing frequency bands were eventually converted 
to a base-band signal of 16 MHz band- width, which was sampled 
at the Nyquist rate. The 16 MHz observing band- width was di- 
vided into 256 spectral channels in the FX correlator. In the GMRT 
array combiner (GAG), the signals from selected antennas were 
added coherently to get the phased array output for each polariza- 
tion which was then fed to the pulsar receiver. The pulsar receiver 
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Figure 1. Spectra with the calculated mean flux densities from relative cal- 
ibration (Table[T). Open circles with error bars represent the flux measure- 
ments from this work and the filled circles represents the flux values avail- 
able in the literature. Spectra flattens below 325 MHz, at frequencies like 
157 and 244 MHz. Known values of the mean flux of PSR B0818~41 at 
different frequencies are: 65 mjy at 400 MHz (Hobbs et al 2004), 20 mJy 
at 660 MHz (Qiao et al. 1995), 5.2 mJy and 9 mJy at 1.4 GHz estimated 
by Taylor et al. (1993) and Qiao et al. (1995) respectively. The spectra from 
our work seems to disagree with mean flux of 9 mJy at 1 .4 GHz (Qiao et al. 
1995) and agrees with mean flux of 5.2 mJy at 1.4 GHz (Taylor et al. 1993). 

computes both the auto and cross polarized powers for each of the 
256 channels. These data were further integrated in time and finally 
recorded to disk at a sampling ra te of 0.512 ms. A suitab l e cali - 
bration procedure as described by iMitra. Gupta & Kudale 1 1 I2OO5I) 
was applied to the recorded data to recover the stokes parameters 
I, Q, U, V. During the off-line analysis the raw data were further 
integrated to achieve the final resolution of 2.048 ms. Table[T]sum- 
marises the observations. 

Average pulse profiles at different observing frequencies were 
obtained by dedispersion with a DM of 1 1 3 .4 pc/ cm^ . Power line 
(50 Hz) and its harmonics were the major contributors to the ra- 
dio frequency interferences present in the data. To remove power 
line related interferences, the dedispersed data were put through a 
radio frequency interference filtering routine which detected most 
(but probably not all) of the power line interferences and replaced 
them by appropriate random noise. The dedispersed, interference 
free data stretch were synchronously folded with the topocentric 
pulsar period. 



3 FLUX AT INDIVIDUAL FREQUENCIES AND 
SPECTRUM 

We determine the mean flux of PSR BOS 18-41 for all the observ- 
ing frequencies at all the observing epochs. For some of the epochs 
we had observed flux calibrator sources and for these cases we de- 
termine the flux from absolute calibration. The calibrator fluxes at 
90 and 20 cm available in the VLA catalog are used as the calibra- 
tor fluxes at 325 and 1060 MHz respectively. However, the calibra- 
tor flux at 610 MHz is obtained by linear interpolation between 90 
and 20 cm fluxes from the VLA catalog and the calibrator fluxes 
at 157 and 244 MHz are obtained by linear extrapolation of the 90 
and 20 cm fluxes. For the epochs where we do not have any flux 
calibrator source observed, we follow a relative calibration pro- 
cedure, using the knowledge of the observing parameters and the 



background sky temperature near th e source (similar to as done in 
iBhattacharvva Gupta & Giil ( I2OO8I) ). As a crosscheck we compare 
the fluxes where both absolute and relative calibrations were pos- 
sible. The flux values determined from absolute and relative cali- 
bration are similar for all the frequencies, except at 157 and 244 
MHz. The reason for the difference in flux values between abso- 
lute and relative calibration could be due to the error in the esti- 
mation of calibrator flux or inaccurate knowledge of the sky tem- 
perature at the position of the pulsar, at the lower frequencies. In 
Table [T] we list the mean flux, for all the five frequencies of obser- 
vations. Though the scintillation properties of PSR B0818— 41 are 
not studied, still we try to understand the effect of scintillation on 
the observed flux, by comparing with two nearby pulsars for which 
the scintillation properties are known. The mean flux densities for 
two nearby pulsars, B0628-28 and PSR B0834-I-06, a re observed 
to ch ange by ~ 30% between the observing epochs jBhat et alTI 
I1999I) . Between different observing sessions of PSR B0818-41, 
we observe 15 to 20% change of the mean flux, for a particular 
frequency (Table|TJ. Hence, it is likely that the estimated mean flux 
values of PSR BOS 18— 41 are affected by the scintillation to a mod- 
erate extent. Fig.[T]presents the spectral behaviour of the mean flux 
densities obtained from relative calibration. The error in mean flux 
is calculated from the standard deviation of the available flux val- 
ues at a given frequency. Our estimation of the mean flux values 
of PSR B0818— 41 at different frequencies are in the same ball- 
park as the earlier results available at 400, 660 and 1400 MHz (Fig. 
[T](. Monotonic increase in pulsar flux is observed with decreasing 
frequency starting from 1400 to 325 MHz with average value of 
spectral index a ~ —2 (where a is defined as / oc 1/"), which is 
similar to the mean value of the p ower law index for a typical pul- 
sar spectrum jMaron et al. II2OO0I) . However, we note that there is a 
tendency of flattening of the pulsar spectrum at frequencies lower 
than 325 MHz (at 244 or 157 MHz). We believe that this flattening 
of the spectrum is associated with the turn-over at low frequencies. 
But because of the uncertainty in the observed flux values, it is dif- 
ficult to estimate the exact frequency of the turn-over. Usually pul- 
sars exhibit such a po wer law spectrum with a low frequency turn- 
over near ~ 1 00 MHz l lMalofeev et al. Il994l) . For PSR BOS 1 8 -4 1 , 
the flux of leading outer, trailing outer and inner region evolve dif- 
ferently with frequency (Fig. Relative strength of the leading 
peak as compared to the trailing peak increases with frequency. 
Inner saddle gets more filled with increasing frequency which is 
evident from the average profiles plotted in Fig. (2] Hence, spec- 
tral indices of different profile components are different for PSR 
BOSIS- 41. Component dependent frequency evolution of flux is 
obs erved for many o ther pulsars (e.g. for PSR Bl 133+16 reported 
by l lBhat et al. lim^ ). 



4 AVERAGE PROFILE EVOLUTION WITH 
FREQUENCY 

Fig. [2] presents the plot of the average profiles at five different ob- 
serving frequencies. We observe a steep sided double peaked av- 
erage profile at 325 MHz, as reported in Paper 1. At 157 MHz, 
the trailing peak is broadened, probably because of scattering ef- 
fects, which become more prominent at lower frequencies. The two 
peaks of the average profile, at 157 MHz, are of unequal intensity, 
leading peak having about half the intensity of the trailing peak. 
With increase in frequency, the leading peak becomes more intense 
and the saddle region joining the two peaks becomes more filled 
up. At 1060 MHz, both the peaks of the double peaked profile are 
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Figure 2. Frequency evolution of the average profile of PSR B0818— 41. Different panels plot the average profiles obtained from the GMRT observations at 
different observing frequencies, 157, 244, 325, 610 and 1060 MHz respectively. Pulse amplitude is in mJy. The separation between two peaks of the double 
peaked pulse profile decreases with increasing frequency following the radius to frequency mapping. 



of comparable intensity. Separation between the two peaks of the 
double peaked average profile (A*!") reduces with increasing fre- 
quency (Table|2]l. A"!? varies with frequency as A<I> oc v'^, with 
average p value of 0.2 (Table |2}. This follows the general trend 
of the "radius to frequency mapping", seen in a majority of pul- 
sars. For another wide profile pulsar, PSR B0826— 34, we reported 
simil ar frequency evolution of A"I> ( iBhattacharvva. Gupta & GiTI 
I2OO8I) . Although the simulations described in Sect. 16.1 1 with both 
the geometries (G-1 and G-2) broadly reproduce the evolution of 
A$ with frequency, it is hard to reproduce the observed intensity 
difference between the leading and trailing peaks. 



5 POLARIZATION STUDY 



Fig.[3]|4]and|5]present the polarization profiles at 325, 610 and 1060 
MHz. Polarization profiles at 610 and 1060 MHz have better signal 
to noise than the published polarization profiles at 660 and 1440 
MHz by loiao et al. and reproduce all the basic polariza- 

tion features reported by them. We have observed PSR B0818— 41 
at several epochs. Analysis of the data from all the epochs produces 
self consistent polarization profiles. The polarization profiles with 
best signal to noise ratio available for a given frequency are pre- 
sented in this work. 
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Table 2. Frequency dependence of A*, P™, Pj" and A^s 
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86.6±0.3 
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.6±1.5Pi 


16.5 






4.3 


10.8 


14.1 


1060 


63±0.3 


64 
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a : simulation is done with G-1 (a =ll°and /3 = -5.4°), with emission height formula rem = 41i?j/~°-^''P_i5° °^^Pi'^' 
b : simulation is done with G-2 (a =175.4°and /3 = -6.9°), with emission height formula rem = 41Ri/~° l*P_i5° °^Pl 
c : error associated in determination of A<l>s is ±1.3 ° 



5.1 Linear polarization 

Significant amount of linear polarization is observed at all the fre- 
quencies (about 45% at 325 MHz, 32% at 610 MHz and 20% at 
1060 MHz) : the percentage of linear polarization is more at 325 
than at 610 and 1060 MHz. This agrees wel l with the percentage 
of linear polarization quoted by lOiao et aT ] ( Il995h . The average 
profile for linear polarization follows the general double peaked 
structure of the total intensity profile. But the separations between 
the two peaks of the linearly polarized profile are, 79° at 325 MHz, 
7 Tat 610 MHz and 60° at 1060 MHz, which are less than the sep- 
arations between the total intensity peaks at individual frequencies. 
We note that this phenomenon appears to have a frequency depen- 
dence: width of linearly polarized and the total intensity profile dif- 
fering by ~10% at 325 MHz whereas by ~5% at 610 and 1060 
MHz. At 325 MHz the leading and the trailing peaks of the average 
linear polarization profile are of similar intensity, whereas in the 
total intensity profile the trailing peak is significantly more intense 
than that of the leading one. At 610 and 1060 MHz the leading and 
the trailing peaks of the linear polarization profile are of similar in- 
tensity, and the same is true for the total intensity profiles. Average 
linear polarization falls off much faster than the total intensity and 
decreases to zero near the outer edges of the profiles. This effect 
of edge depolarization is less in 1060 MHz as compared to 325 
and 610 MHz. One of the striking features seen in Fig. [3] E] and 
[5] is that the PA curve evolves significantly with frequency. In or- 
der to do a closer comparison, Fig. |6] shows the three PA curves 
in one plot. The PA curve shows linear monotonic variations (with 
some amount of iiTegularities) at 325 MHz. Total PA swing is about 
69°at 325 and about 53°at 610 MHz. The PA curve at 610 MHz is 
quite different than that of 325 MHz. For the left half of the PA 
curve at 610 MHz, nearly linear PA variation with pulse longitude 
is observed, which is similar to 325 MHz. A different slope of the 
PA curve is seen at 610 MHz at the middle of the pulse profile. 
The other part of the PA curve at 610 MHz exhibits kinky structure 
and some signature of reversal of the slope of the PA curve is ob- 
served at the edge of the profile. It i s difficult to explain the total PA 
sweep at 610 MHz with the RVM l lRadhakrishnan & Cooklll969h , 
which is indicative of non dipolar emission geometry. Though the 
PA curve is quite noisy at 1060 MHz, we see a behaviour interme- 
diate to that at 325 and 610 MHz (Fig.[6ll. 

Evidence of orthogonal polarization mode jump is observed 
at the two edges of the profile at 325 MHz (Fig.[3](. While investi- 
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Figure 3. Polarization profile of PSRB0818-41 at 325 MHz. The top panel 
shows the average profile in total intensity (solid line), linear polarization 
(dashed line) and circular polarization (dotted line). The lower panel shows 
the variation of PA as a function of longitude. The amplitude is in arbitrary 
units. 



gating the single pulse PA using scatter plot, we observe orthogo- 
nal polarization mode jump near two edges of the pulse profile at 
both 325 and 610 MHz. In 1060 MHz the signal to noise was not 
enough to observe single pulse PA. Orthogonal polarization mode 
jump observed at the edge of the profile (at 325 and 610 MHz), can 
be the reason of the abrupt decrease of the linear polarization near 
the edge of profile. Presence of the secondary polarization mode 
at the edge of the profile can cause the observed reduced width of 
the linearly polarized profile compared to the total intensity profile. 
Stronger orthogonal polarization mode jump observed at 325 MHz 
than 610 or 1060 MHz, can explain the fact that observed differ- 
ences between the separations between the peaks of linearly polar- 
ized and total intensity profile are more at 325 MHz. For another 
wide profile pulsar PSR B0826— 34 we observed similar orthogo- 
nal polarization mode change at the edge of the profile , resulting in 
sharp fall of the linear polarization teiggs et al. II 19851) . 
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Figure 4. Same as Fig.[3] but at 610 MHz 



quencies, whereas sense of circular polarizati on changes from n eg- 
ative to positive direction at high frequencies jOiao et alTIl 19951) . 

The feature of increased circular polarization near the trailing 
peak at 325 MHz, which is not observed at 610 and 1060 MHz, may 
contribute towards the more intense trailing peak than the leading 
one observed at 325 MHz. 

Sense reversals of the circular polarization are commonly as- 
sociated with the core components. However, for PSR B0818— 41, 
at 1060 MHz, it is associated with the conal components (Paper 1). 
There are some conal double pulsars such as PSR B0826— 34, PSR 
B2048— 72 which show a central reve rs al of circular PA b ut no core 
component at all teiggs et alT] ( Il985l) . iHan et al. I fT99^ ). In such 
cases the sense reversal is associated with con al components rather 
than the core components teiggs etal.lll985l) . 

For PSR B0818— 41 we observe a correlation between the sign 
change of circular polarization and the sense of PA sweep: with 
right h and (negative) circ ular polarization accompanying increas- 
ing PA. lHanetal~l l ll998l) reported similar correlation between the 
sense of PA sweep and the sense of circular polarization, using the 
polarization data of 20 conal-double pulsars. 
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Figure 5. Same as Fig.[3] but at 1060 MHz 



5.3 RVM fit to PA sweep : 

RVM jRadhakrishnan & Cook I [19691) is fitted to the PA curve of 
PSR B0818-41 at 325 MHz. While fittin g we used symmetry 
abo ut magnetic axis, i.e. < a < 180°(e.g |Ev"erett et al. 1 120011) 
and ljohnston et all l l2006h ), rather th an symmetry about the mag - 



netic equator, i.e. < a < 90 (e.g. iLvne & Manchester 



1990l) . 



Rankin 



( 11988 
h993h 



iNaravan & Vivekanand I ( Il982l) . iRankinI 
etc). Our attempt to fit RVM to the PA curve at 325 MHz gives 
best fit values of a ~ 175.4° and /3 ~ —6.9°, which is an outer 
LOS geometry. This means that the spin axis is anti aligned to the 
magnetic axis. The angle between the spin axis and the magnetic 
axis is, 180— 175.4~5.6°, and the LOS makes a circle of radius 
5.6 + /3=12°around the rotation axis. Fig. |7] presents a plot of the 
observed PA sweep at 325 and 610 MHz and the RVM curves for 
two possible geometries : a ~ 11°, /3 ~ —5.4° (the geometry de- 
scribed in Paper 1 which incidentally reproduces the middle part of 
PA sweep at 610 MHz) and q ~ 175.4°, /3 ~ -6. 9° (the geometry 
derived from the RVM fit to 325 MHz PA sweep). 



5.2 Circular polarization 

Though circular polarization observed for PSR B0818— 41 at 325 
MHz is quite less, under the trailing component of the pulse profile 
there is a feature of increased circular polarization. Small amount 
of circular polarization is observed at 610 and 1060 MHz. Circular 
polarization changes sign at 1 060 MHz near the mid dle of the pulse 
profile, which is also reported bv lOiao et al.1 d 19951) at 1440 MHz. 
We note that changing sign of circular polarization is no t observed 
at 325, 610 (this work) and 660 MHz dOiao et al. II 19951) . 

There are two possible origins of circular polarization of 
pulsars: either intrinsic to the emission properties and depen - 
dent on the emission mechan is m (e.^. iMichel et al. ] ( fl987h . 
iRadhakrishnan & Rankin I ( Il99(]|). lOangadhara! d 19971)) or gener- 



ated by propagation effects (e.g. lCheng & Ruderman dl979l) ). Cir 



cular polarization of some pulsars changes with frequency and the 
variation of degree of circular polarization with freque ncy is very 
differe nt from pulsar to pulsar (e.g. refer to Table 5 of iHan et al. I 
PSR J2053— 7200 is an extreme example where sense re- 
versal of circular polarization from the positive to negative direc- 
tion is observed near the intersection of two components, at low fre- 



6 UNDERSTANDING THE EMISSION GEOMETRY OF 
PSR B0818-41 



The possible choices of the emission geometries are constrained by 
the following observed characteristics: 

(i) Observed profile width at individual frequencies (Table |2} as 
well as the profile evolution with frequency (shown in Fig.O. 

(ii) Overall pattern of drifting with simultaneous multiple drift re- 
gions (e.g. Fig.|13b, P™, P'P values and the relative intensities of 
the drift regions. 

(iii) Observed PA curves. 

As described in Paper 1 (see Sect. [TJ, the observed drift 
pattern for PSR B0818— 41 can be interpreted as being created by 
the intersection of our LOS with two conal rings on the polar cap 
of a fairly aligned rotator (a ~11°), with an inner LOS geometry 
(P ~ —5.4°). This geometry was developed to give best fit to 
profile shape and drift pattern at 325 MHz. In this work we find 
that the PA traverse at 325 MHz can not be explained by this 
geometry, but surprisingly the middle part of the PA swing at 610 
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Figure 6. Observed PA curves at 325 (solid), 610 (doted line), 1060 (dash- 
dot) 



100 

80 
60 








40 




'v, 

\. 




20 

-20 
















-40 








-60 








-80 
100 









-60 -40 -20 20 40 60 



Pulse Phase (deg) 

Figure 7. Observed PA curves at 325 (solid circles), 610 MHz (solid trian- 
gles) and the RVM curve for a ~ 175.4°and /3 ~ -6.9°(dotted line) and 
for a ~ 11°, 13 5.4°(dashed line) 

MHz is satisfied by this geometry (Fig. |7}. However, best RVM 
fit to the PA curve at 325 MHz as well as left half of PA sweep at 
610 MHz is obtained with a ~175.4°and P ~ -6.9°(Sect. O. 
Therefore, we have two competing geometries, 
G-1 {a =ll°and P — —5.4°) : geometry used in Paper 1, that 
gives best fit to profile shape and drift pattern at 325 MHz. This 
geometry incidentally produces reasonable fit to the middle part of 
PA sweep at 610 MHz. 

G-2 (a =175.4° and (3 = -6.9°) : geometry obtained by RVM 
fitting to the PA curve at 325 MHz. 

In order to investigate how well these two geometries satisfy 
the above constraints we have exhaustively used the simulation 
tecliniques as described below. 

6.1 Simulation 

The simulations with two rings of emission having same carousel 
rotation period (P4) and the number of sparks (Nsp), with 180° out 
of phase emission between the rings, successfully reproduces 
the PLR in the observed drift pattern. Considering the measured 



=17.5°for the inner region and assuming ~ P2 (true lon- 
gitude separation between adjacent drift bands) for nearly aligned 
rotator, we have, Nsp '-^ 19 — 22. We simulat e the radiation 
from PSR BOS 18-41 using the method described in lOil & Sendvkl 
( I2OOOI) for both G-1 and G-2. While doing the simulation we use 
20 equispaced sparks in each of the two concentric rings, rotating 
with a drift rate of 19.05°/Pi (as ^ I8P1, Sect.[7l (. The em- 
pirical emission altitude formula bv lKiiak & Gil 1 1 (200 3l) is used for 
determining the LOS cut through the beam at different frequencies. 

The gray scale plots of the simulated single pulses for G-1 
and G-2 (Fig[8]and FigllOt reproduces the overall pattern of drift- 
ing. Radii of the inner and outer rings in Fig[8](simulated with G-1) 
are 0.6 Vp and 0.8 Vp, which are identical to the values used in Paper 
1. However, for Fig. [TO] (simulated with G-2) radii of the inner and 
outer rings are 0.7 rp and 0.9 r-p respectively. Ratio of spark inten- 
sities of the inner and outer rings, required to reproduce the relative 
intensities of the inner saddle region and the trailing peak, is ~ 3. 
This is identical to the relative intensities of the inner saddle region 
and the trailing peak of the observed pulse profile at 325 MHz (Fig. 
O. From the simulated single pulses with G-1, we calculate P2" ~ 
16° for both inner and outer regions, which is close to the observed 
Pa™ in the inner region (Pa'" ~ 17.5° in SecO whereas for G-2, 
P2" ^ 14° for the inner and outer regions, which is a bit lower than 
the observed P]" for the inner region. Fluctuation spectrum of the 
single pulses generated from both the geometries exhibits a strong 
feature at ~ I8P1, which is identical to the determined P™ value. 
Simulated single pulse drift pattern with G-1 in Fig[8]have single 
drift bands in the outer drift region on both sides of the inner re- 
gion, which matches with observation. But for the simulated single 
pulse drift pattern with G-2, at times, two drift bands in the outer 
region on both sides of the inner region are observed. The simulated 
average profiles for both the geometries and the observed average 
profiles at 325 MHz are over plotted in Fig. |9] In the simulated 
average profiles the inner saddle region get more filled at higher 
frequencies, which reproduces the observations (Fig.|2ll. Though, 
general double peaked structure of the profiles, profile width at in- 
dividual frequencies as well as the profile evolution with frequency 
are well reproduced by the simulations with both the geometries 
(Table|2]l, features like differences in amplitude of the leading and 
trailing peaks are not achieved with simulations. 

To summarise, we have simulated the radiation pattern of the 
pulsar for both G-1 and G-2. Both the geometry satisfy constraints 
(i), (ii) and (iii) discussed in Sect.|6] Though the overall PA curve 
at 325 MHz, and left half of the PA curve at 610 MHz is fitted 
with G-2, only the middle part of the PA curve at 610 MHz is fitted 
with G-1. On the other hand, G-1 appears to give a better match 
to the observed P™ values in the inner region. We will refer to the 
simulation results at different stages of this paper. 



7 SUBPULSE DRIFTING AT MULTIPLE FREQUENCIES 



Fig.[T2]shows an example of single pulse data of 200 pulses at 325 
MHz. Individual pulses are seen over a wide longitude in the pulse 
window. Organized drifting near the two edges of the pulse profile 
is observed, in the form of linear drift bands going from the later 
to the earlier longitude, as the pulse number increases (we refer to 
this as negative drifting). Fig.[T3]is the gray scale plot of a sequence 
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Figure 8. Simulation of the subpulse drift pattern with simple dipolar ge- 
ometiy for the case: a = 11°, f3 = -5.4°, drift rate D = 19.05°/Pi, 
Nsp = 20. Corresponding P5 = 8.7Pi. 



i 



7 

I 

2 



Leadin: 



Peak ' o^P'^'''' 




Direction of 
Carousel rotation 



Figure 10. Simulation of the subpulse drift pattern with simple dipolar ge- 
ometry for the case: a = 175.4°, /3 = -6.9°, drift rate D = 19.05°/Pi, 
Nsp = 20, Corresponding P5 = 8.7Pi. 




Pulse Phase (deg) 

Figure 9. Observed and simulated pulse profile at 325 MHz with simple 
dipolar geometry for the case: a = 11°, /3 = —5.4°; drift rateD = 
19.05° /Pi. 
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Figure 11. Observed and simulated pulse profile at 325 MHz with simple 
dipolar geometry for the case: a = 175.4°, /3 = —6.9°; drift rateD = 
19.05°/Pi. 



of 200 pulses showing regular drifting, at 325 MHz. As reported in 
Paper 1, three different drift regions with different drift rates — an 
inner region with steeper apparent drift rate flanked on each side 
by a region of slower apparent drift rate. Multiple drift bands (typ- 
ically 3 to 4) can be seen in the inner region for any single pulse, 
whereas the outer regions show only a single drift band. Further- 
more, the two drift regions maintain a constant phase relationship. 
Though, simultaneous multiple drift regions are more prominent 
at 325 MHz, drifting is also observed at 150, 244 and 610 MHz. 
Fig[T5]shows the gray scale plot of the single pulses at 610 MHz. 
Though the single pulses become quite weak at 610 MHz, drift- 
ing is clearly visible in the gray scale plot both in the inner and the 
outer region. In the inner region, 3 to 4 simultaneous drift bands are 
visible and faint single drift bands can be identified on both sides 
of the inner region. Fig. [14] is the gray scale plot of single pulses 



at 244 MHz. Drifting is observed in the leading and trailing outer 
regions, but no signature of drifting is present in the inner saddle 
region. Even after integrating the data up to a resolution of 8 ms to 
detect possible weak drift bands in the inner drift region, we do not 
see any evidence of that. At 150 MHz faint drift bands are observed 
in the trailing outer region, but no evidence of drifting is observed 
either in the inner region or in the leading outer region. 

Occasionally, either leading or trailing outer region partially 
or totally disappears, while normal emission is observed from the 
other regions. Fading of the leading outer region is more commonly 
observed, e.g. around pulse # 1655 to 1675 in Fig.[T3j 1420 to 1440 
in Fig.[T4]and 470 to 490 in Fig.fTSi 

We also observe frequent nulling in this pulsar. Fig.ll6lshows 
some sequences of pulsar nulls. Duration of the nulls varies from 
a few pulses to few hundreds of pulses. Nulling property of this 
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Figure 12. Single pulse sequence of 200 pulses (pulse # 200 to 400) of 
PSR BOS 18-41 from the GMRT observations at 325 MHz. Pulse emission 
is present for a wide longitude range (~ 150 °). Diift bands in the outer 
leading and trailing side are observed from the stack of single pulses. Mod- 
ulations in the inner region are weaker and are indistinguishable in this plot 
in contrast to the outer regions. 



pulsar will be discussed in detail in a follow up paper. 
7.1 Frequent change of drift rates : 

Though the regular drift patterns observed in Fig.[T3]are quite com- 
mon for this pulsar, we observe frequent changes of the drift rates. 
Some extreme examples of changing drift rates such as transitions 
from negative to stationary or stationary to negative drifting are 
also observed. Because of the observability of simultaneous multi- 
ple drift regions with good signal to noise, changing drift rates are 
more prominently seen at 325 MHz. Due to the absence of drifting 
in the inner region at 244 MHz and weaker drift bands at 610 MHz, 
such regions are harder to identify at these frequencies. We inves- 
tigate changes in drift rates for about 10,000 pulses from two dif- 
ferent epochs of observations at 325 MHz. In addition to frequent 
occurrences of small changes in the drift rate, we could find, seven 
transitions from the negative to stationary drift rates, five transitions 
from stationary to negative drift rates, and two possible signatures 
of curved drift bands. We also observe one pulse sequence showing 
stationary drifting interleaved between nulling regions. Fig. |16|[77l 
andllSlillustrate some of these phenomenon. Fig.[T6]shows the sin- 
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Figure 13. Gray scale plot of single pulse data of PSR B0818-41 at 325 
MHz, with the average profile shown on top. Fading of the outer left region 
near pulse number range 1655 to 1675, corresponding fading in the outer 
right region might be near pulse # 1670. 



gle pulse gray scale plot for pulse # 2800 to 3200. Almost station- 
ary drift is observed around pulse # 2870 to 2920 after which the 
regular negative drifting starts. Irregular drifting is observed again 
around pulse # 3107. Fig. [17] presents the single pulse gray scale 
plot for pulse # 3600 to 4000. Stationary drift bands are observed 
from pulse # 3665 to 3698 after which regular negative drifting 
starts from pulse # 3702. The apparent drift rate again slows down 
around pulse # 3900, stationary drifting is observed in following 
few pulses and then the pulsar starts nulling around 3913. Fig. 1181 
presents the single pulse gray scale plot for pulse # 4500 to 4900. 
Transition from negative to stationary drift rate is observed around 
pulse # 4545. Stationary drifting is observed from pulse # 4545 
to 4595, after which the pulsar starts nulling and becomes active 
again at pulse # 4603. Fig.[T9]presents the single pulse gray scale 
plot for pulse # 5600 to 6000. Indication of curved drift bands is 
seen around pulse # 5660 to 5720. For the sequences showing tran- 
sitions from negative to stationary drift, nulling starts immediately 
after the stationary drifting. Similarly for the sequences showing 
transitions of stationary to negative drift, the stationary drifting fol- 
lows just after nulling. Almost no situation is seen (apart from some 
indication in curved drift bands) where the drift direction is posi- 
tive. Possible negative-stationary-negative or negative-stationary- 
positive transitions are generally interrupted by nulls. 

To summarise, frequent changes of drift rates — including 
transitions from negative to stationary and stationary to negative 
drift rates — are observed for PSR B0818— 41. We observe transi- 
tions from negative to stationary drift or stationary to the negative 
drift in different pulse sequences at 325 MHz. Such transitions of 
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Figure 14. Gray scale plot of single pulse data of PSR B0818-41 at 244 
MHz. Drifting is observed in the leading outer and trailing outer regions, 
but not in the inner region. Occasional fading of drift bands in the outer left 
region is seen. 



negative to stationary drift are always followed by nulling. After 
nulls the drifting is stationary (or irregular) for few pulse and then 
the transition from stationary to negative drift is observed. 

Interpretations w ith aliasing: According to 
iRuderman & Sutherland! ( 1 19751) model, observed stationary 
drifting for unaliased case will imply that E x B — 0, caused by 
E — 0, and hence no radiation. Consequently it will be difficult 
to explain the observed stationary drifting and the transitions from 
negative to stationary drift rates without aliasing. In the following 
we try to explain t he observed drifting with aliasing. 

According to iRuderman] ( Il976h , for the case of unaliased 
drifting one should observe negative drifting for inner LOS geom- 
etry and positive drifting for outer LOS geometry. Although, for 
G-1 the observed negative drift, can be explained either by una- 
liased drifting or by aliased drifting with even alias order (Paper 
1), observed occasional changes of drift rates can be easily under- 
stood with aliased drifting. Observed negative drift is explained if, 
< P3 < Pi, and the transitions from negative to stationary drift 
rates can be explained by slight decrease of the drift rate, so that P3 
is equal to Pi . The intrinsic drift rate for the unaliased case is pos- 
itive for G-2. Thus to reproduce the observed negative drifting for 
G-2, we need to consider aliased drifting so that Pi < P3 < 2Pi 
and the apparent drift is of opposite sign than that of the true drift 
rate i.e. the true positive drift rate will appear as negative drift rate 
to the observer. Occasionally drift rate slightly increases so that P3 
is equal to Pi, and we observe stationary drift bands. Hence for 
G-1, P3 is slightly less than Pi will give rise to aliased negative 
drifting and transition from the negative to the stationary drifting 
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Figure 15. Gray scale plot of single pulse data of PSR B0818-41 at 610 
MHz showing regular drifting. Drift bands ai'e observed in the inner and the 
outer regions. 



means slowing down of the carousel so that P3 — Pi. For G-2, P3 
is slightly more than Pi will give rise to aliased negative drifting 
and transition from the negative to the stationary drifting means 
speeding up of the carousel. Similarly the observed stationary to 
negative drift rate transitions can be interpreted as speeding up of 
carousel for G-1 and slowing down of the carousel for G-2. 

Frequent occurrence of nulling either before or after of sta- 
tionary drifting indicates that somehow nulling is associated with 
the condition P3 = Pi . After transitions of apparent drift rate from 
negative-stationary direction, we generally observe nulling, imply- 
ing that, pulsar stops emitting after carousel slows down for G-1 
and after the carousel speeds up for G-2 to satisfy P3 = Pi. Quite 
often, recovery of pulsar emission after null, starts from station- 
ary drifting (P3 = Pi) and then the regular negative drifting fol- 
lows. Association of stationary drifting and nulling indicates that 
the mechanism responsible for pulsar drifting is closely associated 
with the pulsar emission mechanism. 

Sign reversals of drifting and curve d drift bands are observed 
for PSR B0826-34 l lBiggs et al. Ill985l) . unlike PSR B0818-41 
where we do not observe positive drifting (except for some indi- 
cations of curved drift bands interrupted by nulling, e.g. Fig.|19ll. 
Similar to our interpretations for explaining the obs erved transi- 
tions o f negative to stationary drift rates in B0818— 41. lGupta et al. I 
l l2004b proposed that the observed drift rate for PSR B0826-34 is 
aliased version of the true drift rate which is such that a subpulse 
drifts to the location of the adjacent subpulse (or a multiple thereof) 
in about one pulse period. They showed that small variations in the 
mean drift rate are enough to explain the apparent reversal of drift 
directions seen in the data. 
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Figure 16. Same as Fig[T3]for pulse # 2800 to 3200. Almost stationary drift 
is observed around pulse # 2870 to 2920. Then the regular' negative drifting 
starts. Frequent nulling is observed. 



iGupta et"aLl l l2004t) proposed that the small variations in the 
mean drift rate caused by very small heating and cooling effects of 
the polar cap can explain the apparent reversals of drift directions 
observed for PSR B0826-34. Similarly, the slight change of the 
drift rate required for transition form negative to stationary direc- 
tions of apparent drifting for PSR B0818— 41 can be explained by 
very small hea ting and cooling effects in the surf ace temperature of 
the polar cap jGil, Mehkidze & GeDpertll2003l) . We note that the 
transitions of drift rates are rather abrupt for quite a few sequences. 
Which can be interpreted by means of quantum jump or quite fast 
change in polar cap temperature, causing P3 slightly less than Pi 
(giving rise to observed negative drift) to P3 — Pi (giving rise to 
stationary drifting) for G-1, and vice verse. 



7.2 Determination of P-^ : 

For the determi nation of P-^, w e use the fluctuation spectrum anal- 
ysis technique ( lBacker|[l970l) . Fig. [20] plots the phase resolved 
fluctuation spectrum for a sequence of 200 pulses showing regu- 
lar drifting in Fig[T3] The strong peak at 18.3 Pi is identified as 
P™. The value of P3" varies slightly for different pulse sequences 
with regular drifting, from minimum of about 16.7 Pi to maximum 
of 21.4 Pi. For a given pulse sequence, P^ appears to be the same 
for leading outer region, inner saddle region and trailing outer re- 
gion. This is true even for the pulse sequences showing in'egular 
drifting. 

Fig. I2T] shows the fluctuation spectrum for the sequence of 
800 pulses with regular drifting at 610 MHz shown in Fig.[T5] We 
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Figure 17. Same as Fig[T3]for pulse # 3600 to 4000. Transition trom neg- 
ative to stationary drift rate around pulse # 3665, almost stationary drifting 
from pulse # 3665 to 3698, regular drift starts from pulse number 3702. 
Drift rate again slows down around pulse # 3900, stationary drifting is ob- 
served in following few pulses and then the pulsar starts nulling around 
3913. Pulsar becomes active at around pulse # 3920. 



observe a strong peak corresponding to P3" at 18.6 Pi, which is 
the same for both the inner and the outer regions. 

Fig.|22]shows the fluctuation spectrum at 244 MHz for the se- 
quence of 500 pulses shown in Fig. [14] We observe a strong peak 
corresponding to P™ at 18.6 Pi, for the leading outer and the trail- 
ing outer regions. Even after integrating the data up to 8 ms, the 
fluctuation spectrum analysis does not show any periodicity present 
in the inner region. 

Table|2]lists the P™ values, at the individual frequencies, cal- 
culated for pulse sequences with regular drifting. Within error bars, 
the value of P3" (~ I8P1) is the same at all the frequencies. This 
su pports the general result a rrived at, for a larger sample of pulsars 
bv IWeltevrede et al. I ( |200^ and (2007) at 20 and 92 cm respec- 
tively. 



7.3 Determination of P^ : 

From the gray scale plot of the single pulses at 325 MHz (e.g. 
Fig|13ll, it is evident that the P™ values for the outer regions are 
more than that of the inner region. Also it appears that the drift 
bands in the trailing outer region are somewhat steeper than the 
drift bands in the leading outer region. P™ value in the inner drift 
region at 325 MHz is calculated from the auto correlation func- 
tion of the single pulses averaged over the number of available 
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Figure 18. Same as Fig[T3]for pulse # 4500 to 4900. Transition from neg- 
ative to stationary drift rate is observed around pulse # 4545. Stationary 
drifting is observed from pulse # 4545 to 4595, after which the pulsar starts 
nulling and becomes active again at pulse # 4603. 

pulses. The secondary peak of the correlation function will corre- 
spond to the correlation between adjacent drift bands and will give 
us the P2" value. From the location of this secondary peak, 
value for the inner region is found to be 17.5±1.3°for 200 pulses 
(pulse # 1600 to 1800) showing regular drifting (Fig|13t. We de- 
termine similar Pj" values, when calculated for another pulse se- 
quence of 200 pulses showing regular drifting as well as for the full 
data set containing all the 6600 pulses. Similar correlation anal- 
ysis performed for the inner region at 244 MHz, even with 8 ms 
integration, does not show any secondary peak. Due to poor sig- 
nal to noise, the secondary peak position at 610 MHz can not be 
clearly determined from zero pulse offset correlation. From the one 
pulse offset correlation (i.e. while computing the auto correlation 
fu nction, we correlated each pulse to the adjacent pulse, as done 
in lBhattacharvva et al. I J2007h ). we determine, P2'"=16.5±1.3°for 
the inner drift region at 610 MHz. 

Since simultaneous multiple drift bands are not visible in the 
outer regions, P™ values can not be estimated using the method 
described above. Instead we use the following technique to deter- 
mine the PJ" value of the outer regions. The correlation function, as 
described above, is calculated for different pulse offsets. The peak 
of the correlation function with zero pulse offset is at zero lag and 
with increase in pulse offset the peak shifts to higher lags in pulse 
longitude (e.g. Fig.|23|l. We plot the position of the peak of the cor- 
relation function versus the pulse offset. This plot is a straight line 
whose slope should reproduce the slope of a drift band (e.g. Fig. 
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Figure 19. Same as Fig[l3]for pulse # 5600 to 6000. There is some indica- 
tion of curved drift band around pulse # 5660 to 5720. 
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Figure 20. The contour plot of the power spectrum of the flux at 325 MHz 
as a function of pulse phase for the sequence of 200 pulses showing regular 
drifting as in Fig |I3l The left panel shows the power spectrum integrated 
over the entire pulse longitude. The upper panel shows the power integrated 
over fluctuation frequency. The peak fluctuation feature is at 18.3 Pi. 

I24t . From the slope of the fitted straight line and the P™ values, 
one can estimate the P™. We refer to this method as pulse offset 
correlation method (POCM hereafter). A pulse sequence showing 
regular drifting (pulse # 1600 to 1800 shown in Fig|13t. is consid- 
ered for determination of P™ values. As a crosscheck, we apply 
POCM to determine the P2" value of the inner region. The Pj" 
values for the inner drift region at 325 and 610 MHz determined 
using POCM are identical to the respective values determined from 
the secondary peaks of the correlation function. Fig. |23]is the plot 
of the auto correlation functions of the leading outer region calcu- 
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Figure 21. Same as Fig|20] but at 610 MHz for the sequence of 500 pulses 
in Fig |15l The peak fluctuation feature is at around 18.6 Pi. 
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Figure 22. Same as Fig|20] but at 244 MHz for the sequence of 500 pulses 
in Fig |14l The peak fluctuation feature is at 18.6 Pi. 

lated for different pulse offsets at 325 MHz. The peak of the cor- 
relation function shifts with increasing pulse offsets. The positions 
of the peaks of the correlation function versus the pulse offset are 
plotted for the leading outer, inner and trailing outer drift regions 
in Fig. |24] From the slopes of the straight lines fitted to those, we 
determine P™=21.6°for the trailing outer region and P™=3I°for 
the leading outer region. Determined P™ values for leading outer, 
inner and trailing outer regions calculated for another sequence of 
200 pulses showing regular drifting, are also in the same ballpark. 
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Figure 23. Auto correlation results for the leading outer region at 325 
MHz. Curves are for 4 different pulse offsets : soHd curve (pulse offset=0), 
long dash (pulse offset=— 1), short dash (pulse offset=— 2), dotted (pulse 
offset=— 3), dash-dot (pulse offset=— 4). 
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Figure 24. Plot of position of the peak of the correlation functions ver- 
sus the con'esponding pulse offsets. Squares and dashed line : for the inner 
region and the corresponding fitted straight line with slope — 0.96°/Pi ; cir- 
cles and dotted line : for the trailing outer region and the corresponding 
fitted straight line with slope — 1.2°/Pi; triangle and dash-dot line : for the 
leading outer region and the corresponding fitted straight line with slope 
-1.7°/Pi 

We apply this method for determination P™ at 244 and 610 MHz, 
and obtain, P2'"=16.9°for the leading outer, P2'"=29.8°for the trail- 
ing outer regions at 244 MHz and P2"=16.5°for the inner region 
at 610 MHz. But surprisingly, when calculated for different pulse 
offsets for the outer regions at 610 MHz, very less shift (for trail- 
ing outer region) or almost no shift (for leading outer region) of the 
peak of the correlation function is observed which may imply that 
the longitudinal extents of the outer drift region is quite small at 
610 MHz. So we were unable to determine P™ values for the outer 
regions at 610 MHz. The P™ values for the leading outer and the 
trailing outer region at different frequencies are listed in Table|2] 

The absence of the inner drift region at 244 MHz can be in- 
terpreted as, the LOS missing the emission from the inner ring, 
which is supported by the results from the simulations (Sect. |6Tt . 
PJ" value, for the inner drift region at 325 MHz, determined in this 
work is close to the value estimated in Paper 1. P™ in the outer 
regions are more than that of the inner region, which is also evi- 
dent from the visual inspection of the drift pattern at 325 MHz. The 
P™ for the leading outer drift region is more than the P™ for the 
trailing outer region, both at 244 and 325 MHz (Table |2j. This is 
difficult to explain with a concept of circular outer ring around the 
magnetic pole with isotropic distribution of sparks. With the sim- 
ulations (Sect. 16. U . we obtain identical P™ values for the leading 
and trailing outer regions as well as for the inner region. Spark pat- 
terns centered around a local pole rather than the magnetic dipole 
axis may provide a possible explanation. 



7.4 Determination of the subpulse width (A$s) 

To determine the subpulse width (A$s), we compute the autocor- 
relation function, with zero pulse offset, and find the full width at 
half maximum (FWHM). Values of A$a, calculated for a pulse 
sequence with 200 pulses showing regular drifting at 325 MHz, 
are equal to 12.3°, 6.3°and 10.8°for the leading outer, inner and 
trailing outer regions respectively. A$s values calculated for an- 
other pulse sequence with 200 pulses showing regular drifting, as 
well as for the full data set containing 6600 pulses at 325 MHz, 
are identical. We calculate the A$a values for the inner and the 
outer regions at 610 MHz and for the outer regions at 244 MHz. 
Table|2]lists the values of A$s for all the three drift regions deter- 
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Figure 25. Gray scale plot of 27 single pulses from PSR B0818-41 at 325 
MHz. The horizontal an'ow is drawn to illustrate that the peak emission 
from trailing outer and leading outer region are not in phase (offset by ~ 9 
Pi). 
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Pulse number 

Figure 26. Pulse energy under the trailing outer region (solid line) and the 
leading outer region (dashed line) for pulse # 250—350 at 325 MHz. 



mined from our observations at different frequencies. For both 325 
and 610 MHz, A$s values in the inner region are less than that 
of the outer regions. A$s values in leading outer region are more 
than that of trailing outer region at 325 and 610 MHz. However, 
at 244 MHz, A"l>s values for the leading outer and trailing outer 
regions are comparable and more than the corresponding values at 
325 MHz. Evolution of A'I's values with pulse longitude follows 
the trend of evolution of . Though within error bars, A<3>s value, 
for the inner drift region at 610 MHz is less than the corresponding 
value at 325 MHz, which is expected by the radius to frequency 
mapping. Hence, evolution of A"l>s and values with frequency 
and longitude mostly follow similar trend. 



8 PHASE RELATION BETWEEN MULTIPLE DRIFT 
REGIONS 

Simultaneous multiple drift regions maintain unique phase relation- 
ship for PSR B0818— 41 : (1) peaks of the emission from outer 
leading and trailing regions are at constant phase offset (observed 
at 244, 325 and 610 MHz) and (2) emission from the inner and 
outer regions are always locked in phase (observed at 325 and 610 
MHz). These phase relationships are clearly observed in Fig.|25] a 
zoomed plot of 27 single pulses at 325 MHz. 
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Pulse offset 

Figure 27. Cross correlation of pulse energy in the leading and traihng outer 
regions as a function of pulse offset, calculated for pulse sequences showing 
regular drifting. Solid line (325 MHz), long-dash (610 MHz), short-dash 
(244 MHz). 

8.1 Investigation of the constant phase offset between the 

peaks of emission from leading and trailing outer regions 

As seen in the single pulse gray scale plots (e.g. Fig.l25t. the peaks 
of emission from the leading and trailing outer regions as a func- 
tion of pulse number are not in phase and are offset by a constant 
pulse period, which we denote as, Ps^ 9Pi from eye estimation. 
We use the following procedure to study this phase relation more 
quantitatively. For each pulse we compute the total energy under 
the trailing and the leading peaks (up to a point where the energy 
goes down to 50% of the peak value on both sides of the peak). Fig. 
1261 shows the variation of the energy of the leading and the trail- 
ing peaks with pulse number, for a sequence of 100 pulses showing 
regular drifting. Peaks of the emission for leading and trailing outer 
regions are observed to be offset by P5 ~ 9Pi . To get a more ac- 
curate estimation of P5, we cross correlate the energy under the 
trailing and leading peaks for different pulse lags. A similar exer- 
cise is performed at 244 and 610 MHz. Fig. 1271 presents the plot 
of the correlation between the energy of trailing and leading outer 
regions as a function of pulse lag, for 244, 325 and 610 MHz. Ini- 
tially, the emissions from leading and trailing peaks are anticorre- 
lated; the correlation builds up with increasing pulse offsets and 
reaches to maximum value near ~ 9Pi , after which the secondary 
peaks are observed at I8P1 interval. To avoid contributions from 
frequent nulls and irregular drifting, this exercise was performed 
on sequences of pulses showing regular drifting (e.g. Fig.|13t. We 
repeat this exercise for 13 such pulse sequences, each containing 
about 200 pulses showing regular drifting. P5 varies slightly, from 
minimum of 6P1 to maximum of lOPi, with an average of 9.2Pi. 
No apparent frequency dependence of P5 is observed. 

Solution of the aliasing problem : Drifting subpulses are an ex- 
cellent example of under sampling. The aliasing problem is solved 
for only a few drifting pulsars (e.g. PSR B0809+74, B0943+10 
etc) using extra information from specific features present in the 
data. In the following we describe a scheme for solving the aliasing 
problem using knowledge of P5. 

Fig.|28]illustrates the configuration of sparks in two concentric 
polar rings and the possible LOS cuts, for both G-1 and G-2. Here 
Aa is the profile width measured with respect to magnetic axis and 
As = 360/A'^ap, is the angular separation between two consecutive 
sparks measured around the magnetic axis. Now we consider the 
ratio. 
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Figure 28. Arrangements of sparks in two rings in the polar cap and possi- 
ble LOS cuts for G-1 (left panel) and G-2 (right panel). 



Act 



(1) 



where Ii and /i denote the integer and fractional parts of the ratio. 
This gives the total number of sparks occupying the outer ring be- 
tween the points of intersection with the leading and trailing peaks 
of the profile. Here, /i is a measure of observed phase offset : 
/i = will give drift bands which are in phase, whereas /i = 0.5 
will produce a phase offset of 9Pi. The above analysis assumes 
that the spark pattern drifts by a negligible amount in the dura- 
tion of the LOS traverse from the leading to trailing peak. In real- 
ity, the spark pattern will drift by _D/360° during this interval, 
where D is the amount by which the sparks drift in one pulse pe- 
riod (P3" = 360/D) and A<j!> is the separation between two peaks 
of pulse profile. 



Act 



+ 



A.^ D 
360 



(2) 



As As 

For a finite /2, the peak emission from leading and trailing outer 
regions will be out of phase by a constant pulse offset (P5), which 
will be determined by the time taken by the sparks to traverse (1 — 
/2)As°. 

(l-/2)As 



Ps = 



Da 



(3) 



where. Dap (apparent drift rate) is the amount by which a spark 
pattern with drift rate D will appear to move in one pulse period and 
is given by Dap = mod (fcD- As), k = JiVT(n + 1)/2, n being 
the alias order (Paper 1). With the knowledge of viewing geometry. 
Nap, alias order n and drift rate D, one can easily calculate the 
expected P5 using Eqn. [3j which can then be compared with the 
observed value. Therefore, knowing viewing geometry, Nsp and 
observed value of P5, it is possible to constrain drift rate D and 
possible alias order. This technique can solve the aliasing problem 
and will provide an independent estimation of P3 and P4. The only 
missing term is an estimate for Act. This is obtained from the pulse 
width A(j>, using the known geometry. The general formula that 
relates the azimuthal angle measured with respect to the magnetic 
axis (cr) with the rotational phase ((fi) is, 

sm(a + /3) sin{(j}) 



ain{o) — 



(4) 



sm(r) 

where F is the angle between the magne tic axis and the LOS at 
the rotational phase jGupta et al. Il2004) . Fig. 1291 shows the vari- 
ation of (T as a function of (j) around the magnetic meridian ((j>=Q 
and (T=0) for G-1 (solid curve) and G-2 (dashed curve). For both 



Pulse Phase (deg) 

Figure 29. Variation of u with tf) (see Eqn|4]of text) for G-1 (a = 11°, 
/3 = -5.4°, solid curve) and G-2 {a = 175.4°, /3 = -6.9°. dashed 
curve). 



G-1 and G-2, A$ is equal to 88° at 325 MHz (Table|2ll. This corre- 
sponds to Act ~ 58°for G-1 and Act - 130°for the G-2 (Fig.^. 
We consider two concentric rings each with 19 sparks. Using Eqn. 
□ for G-1, we have Act/As = 3.08, i.e. /i = 3 and /i = 0.08. 
Fractional part /i is very small. Therefore, for negligible drifting, 
as the LOS passes from the first to the third spark, it traverses from 
the leading to the trailing peaks. In this case the peak emission 
from the leading and trailing outer regions will be in phase, which 
is also reproduced with simulations with G-1. Similarly for G-2, 
Act/ As — 6.86 i.e. Ii — 6 and /i = 0.86. Therefore, with negli- 
gible drifting, passage of LOS from the first to the sixth pulse and 
an extra amount, /i As°, will imply the traverse from the leading to 
the trailing peak which is also reproduced in simulations with G-2. 

We calculate pulse offset Ps from Eqn. [3] for both the geome- 
tries for the cases of unaliased drifting and aliased drifting with 
different Nsp values ranging from 17 to 22. Table [3] lists differ- 
ent combinations of D, n, Nsp values required to reproduce the 
observed drift pattern and the corresponding estimations of P5 for 
G-1. Table|4]lists the same for G-2. The results presented in Tables 
[3] and |4] illustrate that it is not possible to reproduce the observed 
phase relationship between the leading and the trailing peak with 
unaliased drifting — it is too slow to move the spark pattern so 
much that the emission is correlated after ~ 9Pi . Hence, the drift- 
ing is aliased for PSR B0818— 41. For a particular geometry and a 
choice of Nsp, observed P5 can be reproduced with certain choice 
of D and n. Knowing D and n, one can calculate the P3 using Eqn. 
1 of Paper 1. Since P4 = Nsp P3; P4 is automatically constrained 
for all the cases considered in Table [3] and |4] Cases for which P5 
is close to the observationally estimated value (P5 = 9 ± 2) are 
indicated by bold face in Table |3] and g] For G-1, P5 = 9 ± 2 is 
satisfied for the following choices, 

(1) Nsp = 20, 1st order aUas (n = 2, fc 
P5 = 8.69Pi, P4 ~ I8P1 ~ 10 s 

(2) Nsp — 19, 2nd order alias (n — 4, k 
Ps = 7.4OP1, P4 - 9Pi ~ 5s 

(3) Nsp — 18, 2nd order alias (n — 4, k 
Ps = 10.89Pi, P4 ~ 9Pi - 5 s 

For G-2, P5 = 9 ± 2 is satisfied for the following choices, 

(1) Nsp = 20, 1st order aUas (n = 1, fc = 1), D = 19.05°/Pi, 
Ps = 8.74Pi, P4 - I8P1 ~ 10 s 

(2) Nsp = 17, 2nd order ahas (n = 3, fc = 2), D = 43.40°/Pi, 
Ps = 7.IIP1, P4 ~ 9Pi ~ 5s 



= 1), D = 19.05°/Pi, 
= 2), D = 38.95°/Pi, 
= 2), D = 41.05°/Pi, 
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According to the above analysis, 1st order aliasing with D = 
19.05°/Pi and Nsp = 20 give the best match of P5 with obser- 
vations, for both G-1 and G-2, which is rather coincidental to us. 
Hence, it does not help to distinguish between G-1 and G-2. We 
note that P4 determined in this work is identical to the estimation 
in Paper 1 and supports the interesting result that the carousel rota- 
tion rate is the same as the measured P™. 

We simulate the radiation for each geometry for different com- 
binations of Nsp and D (refer to Sect. |6.1| ). Next we correlate the 
emission under the leading and trailing peaks of the simulated sin- 
gle pulses at different pulse lags (as is done for the real data). P5 de- 
termined from the simulated data is found to be identical to that cal- 
culated from Eqn. [3] (listed in Tables[3]and|4l(. Hence, our method- 
ology for determining P5 is confirmed by simulations. 

To summarise, the observed phase relationship between the 
peak emission from the leading and trailing outer region is an 
unique feature of PSR B0818— 41 and acts as an additional infor- 
mation other than Pi, P2, P3 and P4. Investigation of this unique 
phase relationship can solve the aliasing problem for this pulsar. It 
shows that drifting has to be aliased in PSR BOS 18—41. We provide 
a few combinations of feasible values of the Nsp, D and n that 
can produce the observed out of phase relation. The results from 
the corresponding simulations support our arguments. We propose, 
P4 ~ I8P1 ~ 10 s (for n = 2 for G-1 and n = 1 for G-2). 
The attempt to determine n and P4 has been successful only for 
a few other pulsars : PSR B 0809+74 (P4 = 165 Pi ~ 200 s, 
n = 2, [van Lee uwen et al. I J2003i) ). PSR B0818-13 (P4 = 30 
s. Ijanssen et^ J2004h) PSR B 0826-34 (P4 = 14 Pi ~ 25.9 



2, lOupta et all ll2004h), PSR B0834 +06 (P4 = 14 



18.9 s. lAsgekar&DeshpandeN2005l)), PSR B094 3-H10 



s, n 

Pi ^ _ 

(P4 = 37 Pi ~ 40 s. w = O lDeshpande & R ankin' ('l9 99|)) and 
PSRB 1857-26 (P4 = 147 Pi ~ 89.9 s, n = 0, Mitra & Rankin I 
( I2OO8I) ). Hence, PSR B0818-41 has the fastest known P4, faster 
than PSR B0834+06. 



8.2 Phase locked relation (PLR) between inner and outer 
ring 

From the gray scale plot of the single pulses (e.g. Fig. [T3} it is 
evident that inner and outer drift regions are locked in phase : in- 
ner drift bands start from almost halfway between two drift bands 
in the leading outer region and always remain connected with the 
drift bands from the trailing outer region. This constant phase re- 
lation does not get perturbed by nulling or irregular drifting and is 
maintained even during the sequences of irregular drifting as well 
as after nulling. The PLR implies that the emission from the inner 
and the outer rings are always in fixed phase. Furthermore, nulling 
and changing drift rates are observed to be simultaneous in both the 
rings. Interestingly, a 180° phase shift between the spark locations 
on the inner and the outer rings is needed in the simulation to match 
the PLR of the observed drift pattern of PSR B08 18-41. Such an 
arrangement of sparks corresponds to a maximally packed struc- 
ture on the polar cap. Spark discharges occur in every place where 
the potential dr op is high enough to ignite and d evelop pair produc- 
tion avalanche jRuderman & Sutherland 1 1975h . so sparks populate 
the polar cap as densely as possible. At the same time, such an ar- 
rangement of sparks should be quasi rigid, a property which is also 
suggested by the data. 

PSR B0826— 34 is another wide profile pulsar for 
which pre sence of si m ultaneo us multiple drif t region s are 
observed teiggsetal.l (Il985l). lEsamdin et alTI l l2005l) and 
iBhattacharvva. Gupta & Gil ( I2OO8I) ). For this pulsar the main 



pulse and the inter pulse emission are i nterpreted to be co ming 
from two conc e ntric r ings of emiss ion 1 Gupta et al. 1 2004 ) and 
lEsamdin et al. ] J200l ). Fig. 9 of IBhattacharvva. Gupta & Gil I 
1 20081) presents the gray scale plot of the single pulses from PSR 
B0826— 34 at 610 MHz. This plot illustrates simultaneous drifting 
in the main pulse (~ 5 to 6 drift bands) and inter pulse (~ 3 to 
4 drift bands) region. The drift bands in these regions are locked 
in phase implying th at the emission from t he inner and the outer 
rings are in phase. lEsamdin et al. ] 1 I2OO5I) applied the method 
of phase tracking to the single pulse data and detected in total 
13 drift bands in the pulse window. PLR betwee n the inner and 
outer rings is evident from their results (Fig. 6 of lEsamdin et al. I 
((2005!)). For PSR B0826-34 we observe frequent nulling and 
changes of drift rates which are simultaneous for both the inner 
and outer rings. In addition, IBhattacharvva. Gupta & Gil I 1 I2OO8I) 
reported significant correlation between total energy of the main 
pulse and the inter pulse. They found that with increasing pulse 
lag this correlation follows similar trend as the autocorrelation 
of total energy of the main pulse or the inter pulse (Fig. 7 of 
IBhattacharvva. Gupta & Giil 120081) ). The only other pulsar for 
which simultaneous occurrence of multiple drift regions and phase 
relation between the drift regions are observed is PSR B08 15-1-09. 
But the mirrored drift bands observed for PSR BOS 1 5 + 09 can not 
be explained within the lRuderman & Sutherland model. 

Some pulsars - e.g. PSR BI237+25, B1857-26 - clearly 
shows multiple compo nents in t heir p rofile (mostly five compo- 
nents, denoted by "M"). lRankinl l [l993h have recognised these pul- 
sars as a distinct observational class of objects, where the emission 
is coming from the core and two concentric conal rings. No signa- 
ture of drifting is detected in the inner ring for these pulsars, possi- 
bly because of the fact that the inner ring is too weak to detect any 
single pulse behaviour. For example, in PSR B 1237+25 drifting is 
observed in the outer ring with fluctuati on feature at P3 = 2.8Pi 
jZuzana & Rankin 112003) . In Fig. 6 of IZuzana & RankinI | |2005|) 
we observe a faint fluctuation feature around P3 = 2.8Pi corre- 
sponding to the inner ring. 

To summarise, signatures of drifting from more than one ring 
are observed for only a few pulsars and for all of those pulsars, 
emission from inner and outer rings are locked in phase. No counter 
example is observed. This requires common drift rate in the outer 
and inner rings, implying that emission in the two rings are not in- 
dependent, and that the conditions responsible for drifting are sim- 
ilar in both rings. This demands that the conditions are similar ev- 
erywhere in the magnetosphere, which can in turn put constraints 
on the theoretical models attempting to explain the pulsar emission 
mechanism. Our finding of PLR between the emission from the 
inner and the outer rings can be explained considering the pulsar 
emission as a pan magnetospheric phenomenon. 

According to lGil & Sendvk at any given time, the po- 

lar cap is populated as densely as possible. They consider, charac- 
teristic dimension of a spark and typical distance between sparks 
being equal to polar cap height h. They do not distinguish between 
the radial and the azimuthal directions, which may favor similar 
conditions (electric field and magnetic field) for both the rings of 
emissio n, leading to pa n magnetospheric conditions across the po- 
lar gap. IWright I l l2003h proposed an alternative model for the pul- 
sar magnetosphere which considers the magnetosphere as an inte- 
gral whole, and suggests that the inner and outer cones found in 
integrated profiles reflect the two intersections of the null surface 
with the light cylinder and the co-rotating dead zones respectively. 
They consider the pulsar radio emission phenomenon as a global 
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Table 3. Calculation of P5 for G-1 (the inner geometry) for different spark arrangements 



Drift rate 

C/Pi) 


Number of 
sparks 

\/ ' sp ! 


Azimuthal separation 
between sparks 


Acr/As 


Apparent drift 

{Dap) 


Alias 
order 

(n, k) 


/2 


Pulse 
offset 


0.95 


19 


18.95 


3.09 


0.95 





0.10 


17.86 


20 


19 


18.95 


3.09 


1.05 


2, 1 


0.34 


11.82 


38.95 


19 


18.95 


3.09 


1.05 


4, 2 


0.50 


7.40 


57.89 


19 


18.95 


3.09 


1.05 


6, 3 


0.83 


3.15 


76.84 


19 


18.95 


3.09 


1.05 


8, 4 


1.08 


17.46 


21.05 


18 


20 


2.92 


1.05 


2, 1 


1.18 


15.53 


41.05 


18 


20 


2.92 


1.05 


4, 2 


1.43 


10.89 


22.23 


17 


21.18 


2.76 


1.05 


2, 1 


1.02 


19.72 


43.40 


17 


21.18 


2.76 


1.05 


4, 2 


1.26 


14.81 


19.05 


20 


18 


3.25 


1.05 


2, 1 


0.51 


8.69 


37.05 


20 


18 


3.25 


1.05 


4, 2 


0.75 


4.22 


18.17 


21 


17.14 


3.41 


1.05 


2, 1 


0.67 


5.34 


35.34 


21 


17.14 


3.41 


1.05 


4, 2 


0.92 


1.36 



f: combinations producing, P5 = 9 ± 2 are indicated by bold face. 



pan magnetospheric phenomenon. Our results are consistent with 
these kinds of models. 



9 SUMMARY 

Following are the interesting new results from our investigation of 
subpulse drifting and polarization properties of PSR B0818— 41 at 
multiple frequencies: 

(1) We estimate the mean flux of PSR B0818-41 at 5 different 
frequencies and show that the spectrum flattens at frequencies 
lower than 325 MHz (at 244 or 157 MHz), providing indication of 
a low frequency turn-over. 

(2) We present the average profiles at five different frequencies. 
The evolution of profile width with frequency follows radius to 
frequency mapping. 

(3) Significant linear polarization is observed at 325, 610 and 1060 
MHz. Linear polarization profiles follow similar trend as the total 
intensity profile and get depolarized at both edges of the profile. 
This can be explained by the observed orthogonal polarization 
mode jump at the edges of the profile. Polarization angle sweep 
across the pulse profile evolves remarkably with frequency (at 



325, 610 and 1060 MHz), which is not generally observed in 
other pulsars. Though very less circular polarization without any 
signature of changing handedness is observed at 325 and 610 
MHz, circular polarization changes sign at the middle of the pulse 
profile at 1060 MHz. 

(4) In addition to the remarkable subpulse drift observed at 325 
MHz, we report subpulse drifting at 244, and 610 MHz. At 244 
MHz subpulse drifting is observed only in the leading and trailing 
outer regions and not in the inner region. Though the drift bands 
are weaker, subpulse drifting is observed in both inner and outer 
region at 610 MHz. Frequent changes of drift rates, including 
transitions of drift rates from negative to stationary and vice verse, 
are observed for this pulsar. Such transitions of drift rates appear 
to have some connection with nulling. P™, and A^s are 
determined for the inner and the outer drift regions. Though P™ is 
observed to be the same for the inner and outer drift regions, P™ 
and A$s are different for different drift regions. 

(5) We report that the peaks of the emission from the trailing and 
leading outer regions, as a function of the pulse number, are offset 
by a constant interval, P5~9Pi. We also report a phase locked 
relation (PLR) between the inner and outer drift regions for PSR 
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Table 4. Calculation of P5 for G-2 (the outer geometry) for different spai'k arrangements 



Drift rate 

(°/Pl) 
(D) 


Number of 
sparks 


Azimuthal separation 
between sparks 


Acr/As 


Apparent drift 

(£>ap) 


Alias 
order 
(n, fc) 


h 


Pulse 
offset 

v,^ / 


1.05 


19 


18.95 


6.87 


1.05 





0.88 


2.25 


20 


19 


18.95 


6.87 


1.05 


1, 1 


1.13 


15.71 


38.95 


19 


18.95 


6.87 


1.05 


3, 2 


1.37 


11.32 


57.89 


19 


18.95 


6.87 


1.05 


5, 3 


1.62 


6.92 


76.84 


19 


18.95 


6.87 


1.05 


7, 4 


1.86 


2.52 


21.05 


18 


20 


6.51 


1.05 


1, 1 


0.76 


4.48 


41.05 


18 


20 


6.51 


1.05 


3, 2 


1.00 


18.83 


22.23 


17 


21.18 


6.14 


1.05 


1, 1 


0.40 


12.02 


43.40 


17 


21.18 


6.14 


1.05 


3, 2 


0.65 


7.11 


19.05 


20 


18 


7.23 


1.05 


1, 1 


0.49 


8.74 


37.05 


20 


18 


7.23 


1.05 


3, 2 


0.73 


4.56 


18.17 


21 


17.14 


7.59 


1.05 


1, 1 


0.85 


2.43 


35.34 


21 


17.14 


7.59 


1.05 


3, 2 


1.09 


14.73 



f: combinations producing, P5 = 9 ± 2 are indicated by bold face. 



B0818-41. 

Interpretations of the analysis results, aided by the simula- 
tions, lead to following new outcomes, 

(1) Based on the frequency evolution of the average profile, ob- 
served PA swing and results from subpulse drifting, we converged 
on two possible choices of emission geometry: G-1 (a = ll°and 
/3 = -5.4°) and G-2 (a = 175.4° and /? = -6.9°). Pulsar ra- 
diation pattern simulated with both the geometries reproduces the 
observed features, except for some differences. Rather surprisingly, 
all other observed properties are found to be satisfied for both these 
geometries. 

(2) We explain the observed subpulse drifting and transitions of ap- 
parent drift rates from negative to stationary as being due to alias- 
ing. The frequent transitions from negative to stationary drift rates 
observed for PSR B0818— 41 can be understood by slight changes 
in polar cap temperature leading to change in the drift rate so that 
P\ = Pi is satisfied and longitude stationary drifting is observed. 

(3) The absence of the inner drift region at 244 MHz is interpreted 
as, the LOS missing the emission from the inner ring. This is sup- 
ported by the results from the simulations. 

(4) A new technique is introduced in this paper for resolving alias- 



ing, using the constant offset (P5 ^ 9Pi) between the peak emis- 
sion from the leading and trailing outer regions. From the result 
of this technique, we propose that the subpulse drifting for PSR 
B0818— 41 is most likely first order aliased, and the conespond- 
ing carousel rotation period P4 = 10 s. This implies that PSR 
B0818— 41 has the fastest known carousel. 

(5) 180° phase shift between the spark locations on the inner and the 
outer rings is needed in the simulation to match the observed PLR 
of the drift pattern, which implies that the sparks that are maximally 
packed on the polar cap. 

(6) For all pulsars for which we know drifting from more than one 
ring (e.g. PSR B0826-34 and possibly PSR B0815+09), the drift 
pattern in the inner and outer rings are always phase locked. This 
could be a significant constraint for the theoretical models of pulsar 
radio emission, and favors an emission model that considers pulsar 
radio emission as a pan magnetospheric phenomenon. 

Hence, the results from the observations, simulations and the 
follow up interpretations, enable fair amount of insight towards 
the emission properties of the unique drifting pulsar B0818— 41. 
This study will also have significant contribution to constrain the 
theoretical models explaining the possible emission mechanism of 
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pulsars. 
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